MAD/PH/745 
SNUTP 93-6 
YUMS 93-1 
March 1993 



Associated J/ip + 7 production as a probe of the 
polarized gluon distribution 



M. A. Doncheski 
Department of Physics 
University of Wisconsin USA 
Madison, WI 53706 

and 

C. S. Kim 
Department of Physics 
Yonsei University 
Seoul 120, Korea 



Abstract 

Associated production of J ftp and a 7 has recently been proposed as clean 
probe of the gluon distribution. The same mechanism can be used to probe the 
polarized gluon content of the proton in polarized proton-proton collisions. We 
study J ftp + 7 production at both polarized fixed target and polarized collider 
energies. 



Interest in high energy spin physics has been recently revived with the result from 
(and interpretations thereof) the EMC collaboration [[]J on polarized fi — p scattering. 
Processes in polarized pp collisions (such as achievable at an upgraded Fermilab fixed 
target facility or at a polarized collider M) sensitive to the polarized gluon content of 
the proton, such as jets[§, [|, |j, direct photons]!, ||, [7|, and heavy quark production]^], 
have been discussed. Another intriguing suggestion, due to Cortes and Pire|§, is to 
consider X2(cc) production where the dominant lowest-order subprocess would be gg — > 
X2- The partonic level asymmetries for X2/X0 production have been calculated in the 



context of potential modelsjll] and are large. Low transverse momentum quarkonium 



production in polarized pp collisions using other methods has also been considered 0, [12 
as has high p T ip production [|13[ . 

In all cases of charmonium production, the experimental signal is £ + £~ (£ = e or 
/i) with the lepton-lepton invariant mass giving the J/tp mass, since xj can decay 
radiatively to J/if) + 7, and the J ftp signature is quite clean. As has been noted [O] , 
the question of extracting the gluon distribution is made less clean by the multitude 
of contributing processes, e.g.: 



9 + 9 - 


-> Xo,2 


9 + 9 ~ 


-> Xj + 9 


q + 9 - 


■* Xo,2 + q 


q + q - 


-> Xo,2 + 9 


9 + 9 ~ 


-> J/^P + g 


9 + 9 ~ 


-> 6(-> J/ij + X) + b 


q + q - 


-> 6(-> J/^j + X) + b. 



The simplicity of the Cortes and Pire idea is now gone. A full 0(a? s ) calculation of 
the spin- dependent production of \j is necessary. At low p T , \j production will also 
involve q + g and q + q initial states, while at high p T in addition the 2 — > 2 kinematics 
make the extraction of parton distribution functions less direct. Furthermore, a very 
careful calculation is required because even processes with small cross section can have 
a large effect on the asymmetry. The extraction of Ag(x,Q 2 ) using inclusive J/ip will 
be a challenge. 

Recently, J/ip produced in association with a 7 has been proposed as a clean 
channel to study the gluon distribution at hadron colliders [0. The radiative xj 
decays can produce J/ip at both low and high p T , but the photon produced will be 
soft (E ~ 0(400 MeV)). If we insist that the experimental signature consist of a J/ip 
and 7, with large but equal and opposite p T there is only one production mechanism 

g + g -> J/iJ) + 7. (2) 

Following Ref. |15[| , this mechanism has been proposed in Ref. [[HJ to study the po- 
larized gluon distribution in polarized fixed target experiments; we perform a more 
detailed analysis, including the analysis of this mechanism at the Brookhaven Rela- 
tivistic Heavy Ion Collider (RHIC) at both 50 GeV and 500 GeV center of mass energy 
and at the Superconducting Super Collider (SSC). Polarized proton-proton operation 
is being considered for RHIC, for at least several months data collection, while the 
tunnel design of the SSC has been modified for the possible future inclusion of the 
Siberian Snakes needed for polarized proton-proton mode. Also, we list the full set of 
helicity amplitudes for this process, explicitly stating the Lorentz frame in which the 
J/ijj helicities are given. 
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The full helicity amplitudes for g + g — > J jib + 7 can be calculated following 
the approach of Gastmans, Troost and Wu|P-7||, with the addition of explicit helicity 
polarization vectors for the J/ib. A convenient set of polarization vectors can be 
found in Bohm and Sack ||18|| . These polarization vectors reduce to the usual massive 
vector boson (+,—,0) polarization vectors in the parton center of mass frame, and so, 
although the expressions for the helicity amplitudes have Lorentz invariant form, the 
(+,—,0) only refer to the Jftb helicity in this one particular frame. We find only one 
independent helicity amplitude (M(++, ++), where the '++,++' refer to the helicity 
°f 9i92ilJ/'4 ) respectively), and the remaining 5 non-zero helicity amplitudes can be 
found by crossing and parity symmetries: 

M( ++ , ++ ) = M(-,-)=C SiS ~ m 
M (+-,-+) = M(-+,+-) = C 
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M(-+,-+) = M(+-+-)=C - — ^ — — - (3) 



where C = ^ e q e 9sR{fyMb — Here, M is the J lib mass, s, i and u are the usual 
Mandelstam variables, -R(O) is the radial wavefunction at the origin of the cc in the 
J I ip and a, b are the color indices of the incident gluons. Thus, the (spin and color) 



summed and averaged matrix element squared can be found fl5 

(lQir) 2 aa 2 M\R(0)\ 2 ' 



\M(g + g^J/1> + 7 )\< 
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(4) 



(u - M 2 ) 2 (s - M 2 ) 2 (s-M 2 ) 2 (t-M 2 f 
|i?(0)| 2 can be related to the leptonic width of the J/ib: 

T(J/^^e + e-) = ^|i?(0)| 2 = 4.72keV 

(^(O)! 2 = 0.48 GeV 3 . (5) 



We are interested in the longitudinal spin-spin asymmetry, defined as: 

_ <t( + + )-<t( + -) 

L - *++)+*(+-) (6) 

where cr(++) (<r(H — )) is the cross section for the collision of 2 protons with the same 
(opposite) helicities. This can be calculated in the parton model, 

A LL a = J dxi dx 2 d LL a Ag(xi, Q 2 ) Ag(x 2 , Q 2 ) (7) 



where a is the parton level cross section (related to \M\ 2 given earlier), Ag(x, Q 2 ) is the 

polarized gluon distribution in the proton (= (g + (x, Q 2 ) — g~(x, Q 2 )) where g + (x, Q 2 ) 

(g~(x,Q 2 )) is the distribution for gluons with the same (opposite) helicity as that of 

the proton) and clll is the parton level asymmetry 

£(++)-£(+-) 
ttLL ~ a(++) + a(+-y (8) 

Given the known helicity amplitudes for this process, the parton level asymmetry is 

simply 

. s 2 (s - M 2 ) 2 - t 2 (t - M 2 ) 2 - u 2 {u - M 2 ) 2 
aLL ~ § 2 (s - M 2 ) 2 + t 2 (t - M 2 ) 2 + u 2 (u - M 2 ) 2 ' ( ' 

Measurable quantities of interest are the p T distribution and the joint p T — y\ — y 2 

distribution with yi = y 2 = 0, where 1/1(2) is the rapidity of the 7 {J /if}). In the latter 

case, both partons have the same Bjorken-a; (which is a function of p T only). The 

corresponding asymmetries are given by: 

a(++)-a(+-) 
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Note that A\ L is proportional to [Ag(x(p T ), Q 2 )} 2 - Another interesting theoretical 
concept (though not measurable experimentally) is the average clll, or 'resolving 
power'. It is defined in the following way 

(a LL )(j = J dxi dx 2 d LL a g(x x , Q 2 ) g(x 2 , Q 2 ). (11) 

As we wish to determine if a given experimental scenario can shed light on the size of 
the polarized gluon in the proton, we need, in addition to calculating the asymmetry, 
to estimate the experimental uncertainty in the asymmetry. We will approximate the 
uncertainty by the statistical uncertainty, since ratios of cross sections should be rela- 
tively free of systematic uncertainties. The statistical uncertainty in the measurement 
of an asymmetry is given by 5 A, where 

(12) 

Vn 



and N is the number of events. 

We examine this process in several different experimental settings. First, we con- 
sider an hypothetical fixed target experiment and to be specific, take the proton beam 
energy to be 800 GeV (such as would exist at the upgraded Fermilab fixed target fa- 
cility). In order to estimate the luminosity possible at such an experiment, we must 
make some assumptions. First, the Main Injector at Fermilab can provide ~ 10 14 (un- 
polarized protons) /sec, with a 65% duty cycle ||19|| . We'll assume a one month run, at 



a much reduced proton rate (say, a factor of 100), combined with a small polarized gas 
(H2) jet target (approximately 1 cm long). This will give, we think, a very conserva- 
tive estimate of / Cdt = 50 pb -1 . We place no cuts on the rapidity of the photon or 
J/if>, nor on the p T of the photon or leptons. We find a cross section of approximately 
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200 pb, most of which is at low p T . The resolving power (or average clll) is found to be 
about 28%. We use the polarized distributions of Bourrely, Guillet and Chiappetta [pO[| . 
They provide 2 sets of distributions, one with a large polarized gluon distribution and 
small polarized strange quark distribution (we'll refer to it as the set BGCO) and one 
with a moderately large polarized gluon and moderately large polarized strange quark 
distribution (we'll refer to this set as BGC1). The p T distribution is shown in Figure 
la (in cross section) and in Figure lb (in A 2 LL ). We were also interested the asymme- 
try A\ L , (technically, instead of taking y\ and y<i derivatives, we bin the events in the 
usual way, displaying the contents of the bin with —0.1 < yi,y2 < 0.1). The results 
are shown in Figure 3a (distribution in cross section) and 3b {A\ L vs. p T ). We present 
in Table 1 the total number of events expected (at all p T and yi^ consistent with 
our cuts) as well as the 'resolving power' and asymmetry A\ L and an estimate of the 
statistical uncertainty, SA 1 LL . We also list the number of events in a single p T bin (jp T 
given in the table caption), and A 2 LL and $A 2 LL for that particular p T bin. Finally, we 
present the the number of events in the same p T bin, further restricting the events to 
lie within |y 12 | < 0.1, and the value of A\ L and 5A\ L in the particular p T bin. These 
are representative results. Higher statistics can be obtained by the inclusion of all p T 
bins. 



At this point, we would like to further address the work of Ref. [[TjJ. The large 
asymmetries shown are surprising, and in our opinion not correct. The parton level 
asymmetry, making the following replacements for t and u {i.e. working in the parton 
center of mass frame): 

1 

2 
1 



-~(s - M 2 )(l - cos i 



u = - M 2 )(l + cos9) (13) 
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reduces to 



1 - |[(1 + 6 cos 2 9 + cos 4 9) + 2M1(1 - cos 4 9) + - cos 2 #) 2 ] 
1 + |[(1 + 6 cos 2 9 + cos 4 9) + ^(1 - cos 4 <9) + ^(1 — cos 2 9) 2 } 



Here cos# is measured in the parton center of mass frame. It is obvious that for 
cos# = ±1, (111 is a minimum (actually zero), and so, for any s, the maximum of an 
should be at cos 9 = 0. In this limit, the asymmetry reduces to 

1 i( s + M 2 



a LL (cos9 = 0) = X - ( 2 . (15) 

Two further limiting cases are possible, namely production at threshold (s = M 2 ) 
which gives a LL = | and production at very high energy (s — > oo) which gives 
«ll — g- F° r = y/s = 38.75 GeV (the fixed target energy considered both here 
and in Ref. |16|]), the parton level asymmetry is near it's maximum value. Since 
Ag(x,Q 2 )/g(x,Q 2 ) < 1 generally, the maximum observable asymmetry is bounded 
by the maximum parton level asymmetry. Thus we are unable to understand the 
prediction, in Ref. fl6fl , that the observable asymmetry can be as large as 85%. 

Next, we consider collider experiments at RHIC. RHIC is a high luminosity (C = 
2 x 10 32 cm _2 sec _1 = 6000 pb _1 /yr) collider capable of producing proton on proton 
collisions for center of mass energies between 50 and 500 GeV. A program of polarized 
proton on proton collisions, at full energy and luminosity, is being discussed [21]. We 



will assume a nominal running time of 2 months, at full luminosity, for 50 GeV and 
500 GeV each. In order to be somewhat conservative, we will estimate event num- 
bers based on 300 pb _1 integrated luminosity. We will assume a generic collider type 
detector, and in order to simulate the acceptance we will require the photon and elec- 
trons observed to lie in the rapidity range \y\ < 2 (this simulates the acceptance of 



the proposed STAR detector at RHIC 1 22], level 2 for photons and electrons. We will 
not consider the possibility of the detection of the final state at RHIC). Further- 

more, we will (rather arbitrarily) require the p T of the photon larger than 1 GeV in 
the following discussion. We present our results for the p T distribution in Figure 3a, 
and A\ L in Figures 3b (yfs = 50 GeV) and 3c (y/s = 500 GeV). See Figure 4a for 

dp dyidy 2 VS ' Pt and Fi S ures 4b (V^ = 50 GeV ) and 4c (V** = 500 GeV ) for A ll 
vs. p T . The 'resolving power' increases with energy (actually p T ), even though the 

observed asymmetry decreases. This is simply a consequence of the behavior of the 

polarized gluon distribution. Please refer to Table 1 for some representative results. 

Finally, we consider a collider experiment at the SSC. The luminosity of the SSC 
is C = 10 33 cm _2 sec _1 = 30000 pb~ 1 /yr. We will again assume a running time of 
2 months at full luminosity and energy, and conservatively calculate event numbers 
based on 1500 pb~ x integrated luminosity. We require the photons and leptons to have 
p T > 10 GeV and lie in the range \y\ < 2.5 (these approximate the acceptances of 
the SDC detector |53|). In this case, the resolving power is quite high, (d LL ) = 60%, 
although because of the extremely small-x probed the observed asymmetry A l LL is tiny. 
Similarly, A 2 LL and A\ L are both smaller than 1% for all p T < 125 GeV, while there 
will only be a handful of events at (or beyond) p T ~ 25 GeV, so there is no observable 
asymmetry. Again, see Table 1 for some representative results. 

In conclusion, we have studied the process p + p — > J/ij) + 7 + X in polarized 
proton-proton collisions. We first presented the necessary helicity amplitudes and dis- 
cussed the calculation. Then we studied this process at polarized fixed target and in 
colliders, at polarized RHIC (50 and 500 GeV center of mass energy) and at polarized 
SSC. Our results indicate that a polarized (double spin) fixed target program can be 



very useful in the determination of the polarized gluon distribution. It is unfortunate 
that no such experiment is planned. RHIC (especially at lower energies) is an ex- 
cellent probe of the polarized gluon distribution. Since A\ L is directly proportional 
to [Ag(x(p T ), Q 2 )/ g(x(p T ), Q 2 )} 2 , this distribution provides an easy determination of 
the polarized gluon distribution at various x values. It will prove especially useful to 
measure this distribution at several center of mass energies. Even a measurement of 
A 2 LL can provide much useful information (though it is not clear whether the higher 
statistics involved in this measurement will outweight the cleanliness of the extraction 
of the polarized gluon distribution in a measurement of A 3 LL ). The SSC probes a much 
lower x in this process, and since Ag(x, Q 2 )/g(x, Q 2 ) <C 1 there is no measurable asym- 
metry. However, the 'resolving power' at SSC is still very large, so the smallness of the 
asymmetry is purely a consequence of the small-x behavior of Ag(x,Q 2 ). Polarized 
SSC can still be a useful tool for the study of high energy spin properties of the proton 
by utilizing a subprocess that will probe larger x (e.g. heavy Higgs production). We 
should also point out that we have considered only the color singlet model of heavy 
quarkonium production in this paper. A similar analysis can be performed using lo- 
cal duality, if it is determined at HERA that this mechanism contributes to J/ip + 7 
production |24|]. Some slight modifications will be required, namely the inclusion of 



charm in the proton (this effect should be small) and light qq fusion, and in addition 
the modification of the parton level asymmetries. As a final related comment, we plan 
to study J/ijj + 7 production at HERA using a polarized lepton beam and angular 
distributions of the final leptons to learn something of the polarized gluon distribution 
of the photon. 
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RHIC 


86400 


44.7% 


.4% (0.3%) 


4500 


1.7% (1.5%) 


840 


1.8% (3%) 


500 GeV 






.05% (0.3%) 




.2% (1.5%) 




.3% (3%) 


ssc 


8835 


60.2% 


.005% (1%) 
.0006% (1%) 


3000 


.008% (2%) 
.001% (2%) 


540 


.01% (4%) 
.001% (4%) 



Table 1: Summary of representative predictions for J/ip + 7 production in polarized 
proton-proton interactions. Ntot is the total number of events above some minimum 
p T (= GeV for fixed target, 1 GeV for RHIC and 10 GeV for SSC). (a LL ) is the 
'resolving power' as defined in the text (this is independent of the polarized parton 
distributions). A LL and 5A LL are defined in the text; the upper entry corresponds to 
the large Ag(x,Q 2 ) (set BGC0) and the lower entry corresponds to the moderately 
large Ag(x,Q 2 ) (set BGC1). N Pt is the number of events in the particular p T bin 
(0.5-1.5 GeV for fixed target, 1-2 T GeV for RHIC at 50 GeV, 3-5 GeV for RHIC at 
500 GeV and 10-20 GeV for SSC). 
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Figure Captions 



Figure 1 - p T distribution, vs. p T (la) and A 2 LL vs. p T (lb) for large Ag(x, Q 2 
(solid line) and for moderately large Ag(x, Q 2 ) (dashed line) at fixed target. 



Figure 2 " dp% x dy^ =y^ vs - Pt ^ 2a ) and A * LL vs - Pt ^ 2b ) for large A d(x,Q 2 ) 
(solid line) and moderately large Ag(x, Q 2 ) (dashed line) at fixed target. 

Figure 3 - p T distribution, j^- vs. p T (3a) for RHIC at y/s = 500 GeV (solid line) 
and at y/s = 50 GeV (dot-dashed line), and A 2 LL vs. p T for RHIC at y/s = 50 GeV 
(3b) and at y/s = 500 GeV (3c) for large Ag(x, Q 2 ) (solid line) and for moderately 
large Ag(x,Q 2 ) (dashed line). 

Figure 4 - dp ^ dy J yi=y2=Q vs. p T (4a) for RHIC at ^= 500 GeV (solid line) 
and at y/s = 50 GeV (dot-dashed line), and A\ L vs. p T for RHIC at yfs = 50 GeV 
(4b) and at y/s = 500 GeV (4c) for large Ag(x,Q 2 ) (solid line) and moderately large 
Ag(x,Q 2 ) (dashed line). 
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